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ABSTRACT 

Orbital, spin and astrometric parameters of the millisecond pulsar 
PSR J0621+1002 have been determined through six years of timing observa- 
tions at three radio telescopes. The chief result is a measurement of the rate 
of periastron advance, w = 0?0116 ± 0?0008 yr" 1 . Interpreted as a general rela- 
tivistic effect, this implies the sum of the pulsar mass, mj, and the companion 
mass, m 2 , to be M = m x + m 2 = 2.81 ± 0.30 M 0 . The Keplerian parame- 
ters rule out certain combinations of rri\ and m 2 , as does the non-detection of 
Shapiro delay in the pulse arrival times. These constraints, together with the as- 
sumption that the companion is a white dwarf, lead to the maximum likelihood 
values mj = 1.69^30 M© an d m 2 = 0.98^o.'i2 M© (68% confidence). The other 
major finding is that the pulsar experiences dramatic variability in its disper- 
sion measure (DM), with gradients as steep as 0.013 pccm^yr -1 - A structure 
function analysis of the DM variations uncovers spatial fluctuations in the in- 
terstellar electron density that cannot be fit to a single power law, unlike the 
Kolmogorov turbulent spectrum that has been seen in the direction of other pul- 
sars. Other results from the timing analysis include the first measurements of 
the pulsar’s proper motion, // = 3.5 ± 0.3 mas yr -1 , and of its spin-down rate, 
dP/dt = 4.7 x 10 20 , which, when corrected for kinematic biases and combined 
with the pulse period, P = 28.8 ms, gives a characteristic age of 1.1 x 10 10 yr and 
a surface magnetic field strength of 1.2 x 10 9 G. 

Subject headings: stars: neutron — binaries: general — pulsars: individual 
(PSR J062 1 -h 1 002) — relativity — ISM: structure 
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1. Introduction 

Recent pulsar searches have uncovered a new class of binary pulsars. Most 
millisecond pulsars are in nearly circular binary orbits with low-mass He white 
dwarfs, comprising the class of so-called low-mass binary pulsars (LMBPs). The new 
category of systems, the intermediate-mass binary pulsars (IMBPs), have companion 
stars with masses above 0.45 M 0 . Since helium flash occurs at a core mass of 
0.45 M 0 , the stars must be CO or ONeMg white dwarfs. About eight IMBPs are 
currently known. Observationally, they are distinguished by high mass functions, 
fi > 0.015 M 0 ; by moderately spun-up pulse periods, 10 ms < P < 200 ms; and by orbital 
eccentricities that are small, e < 10“ 2 , but somewhat higher than those of the LMBPs 
(Camilo et al. 1996; Camilo et al. 2001; Edwards & Bailes 2001). 

This paper reports on timing measurements of PSR J0621+1002, an IMBP with a 28.8- 
ms spin period in an 8.3-day orbit. The main goal of our observations was to determine the 
pulsar and companion masses through measurement of post-Newtonian orbital parameters, 
particularly the rate of apsidal motion. Measuring apsidal motion is challenging in white 
dwarf-pulsar binaries because their often extremely small eccentricities— values of 10~ 5 
are typical — make it difficult to determine the angle of periastron, and hence to measure 
periastron advance. A comparatively high eccentricity still only e = 0.0025 made the 
detection of apsidal motion feasible for PSR J0621+1002. 

Mass measurements in white dwarf-pulsar systems can be used to constrain theories 
of binary evolution. The LMBPs and the IMBPs share roughly similar histories. They 
both originate when a giant star transfers mass onto a neutron star companion, resulting 
in a spun-up millisecond pulsar and a white dwarf. The histories of LMBPs and IMBPs 
differ, however, in the details of mass exchange. For LMBPs, it is a stable process that 
occurs when the giant swells past its Roche lobe (Phinney k Kulkarni 1994). For IMBPs, 
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it is thought to be an unstable transfer that operates via common envelope evolution 
(van den Heuvel 1994) or super-Eddington accretion (Taam, King, & Ritter 2000). Other 
pulsar binaries, such as double neutron star systems, evolve in still other ways. One way 
to test binary evolution scenarios is by comparing the masses of neutron stars in different 
classes of binary systems to infer differences in amounts of mass transferred. (For a review 
of pulsar mass measurements, see Thorsett & Chakrabarty 1999.) 

The first year of timing observations of PSR J0621+1002 was discussed in Camilo 
et al. (1996), which reported the pulsar’s Keplerian orbital elements, position, and spin 
period. With five additional years of timing data, we have measured the apsidal motion, 
spin-down rate, and proper motion of the pulsar, and we have derived significant constraints 
on Shapiro delay. We also have found sharp variations in the dispersion measure (DM), 
which we use to analyze turbulence in the interstellar medium (ISM) in the direction of 
PSR J 0621+1002. 


2. Observations 

Radio telescopes at Arecibo, Green Bank, and Jodrell Bank recorded pulse times of 
arrival (TOAs) from PSR J0621+1002 on 529 separate days between 1995 March 18 and 
2001 July 1. Table 1 summarizes the observations. The data comprise (1) three nine-day 
campaigns at Arecibo in 1999 May, 2000 May, and 2001 June, each covering a full pulsar 
orbit at two radio frequencies; (2) a handful of additional Arecibo measurements taken 
monthly between 1997 and 2001; (3) irregularly spaced observations at Jodrell Bank on 387 
days between 1995 and 2001, with an average of five days between epochs; (4) twenty Green 
Bank sessions spaced roughly two months apart, each performed over four consecutive days 
at two frequencies; and (5) four campaigns covering the full orbit at Green Bank, two each 
in 1995 and 1998. 
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At the 305-m Arecibo telescope, the Princeton Mark IV data acquisition system 
(Stairs et al. 2000) collected three or four 29-minute data sets each day at 430 and 
1410 MHz. Local oscillators in phase quadrature mixed a 5-MHz passband (10-MHz at 
1410 MHz) to baseband in both senses of circular polarization. The four resulting signals 
were low-pass filtered, sampled, quantized with 4-bit resolution (2-bit at 1410 MHz) and 
stored on disk or tape. Upon playback, software coherently dedispersed the voltages 
and folded them synchronously at the pulse period over 190-second integrations, yielding 
1024-bin pulse profiles with four Stokes parameters. 

The 140-Foot (43-m) telescope at Green Bank observed the pulsar for 30 to 60 minutes 
a day at 370, 575 or 800 MHz. The “Spectral Processor”, a digital Fourier transform 
spectrometer, divided the signals into 512 spectral channels across a 40-MHz passband in 
each of two polarizations. The spectra were folded synchronously at the pulse period over 
an integration time of 300 seconds, producing pulse profiles with 128 phase bins each. 

The 76-m Lovell telescope at Jodrell Bank carried out a typical observation of 30 
minutes at 410, 606 or 1400 MHz. The signal was dedispersed on-line in each of two 
polarizations using filterbank spectrometers with bandwidths of 64 x 0.125 MHz for the 400 
and 600 MHz data and 32 x 1, 32 x 3, and 64 x 1 MHz for the 1400 MHz data. The detected 
signals were folded synchronously to make a pulse profile. 

In all cases, conventional techniques were used to measure pulse arrival times. Spectral 
data were dedispersed and summed to produce a single total-intensity profile for a given 
integration. Each profile was cross-correlated with a standard template to measure the 
phase offset of the pulse within the profile. Different templates were used for each receiver 
at each telescope, and for each spectrometer at Jodrell Bank. The offset was added to the 
start time and translated to the middle of the integration to yield a TOA. In a further step, 
sets of TOAs from Arecibo and Green Bank were averaged over intervals of 29 minutes 
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(sometimes longer for Green Bank) to make a single effective TOA for the interval. Each 
observatory’s clock was corrected retroactively to the UTC timescale, using data from the 
Global Positioning System (GPS) satellites. 


3. Timing Model 

The TOAs were fit to a model of the pulsar’s orbital, astrometric and spin-down 
behavior using least-squares methods. The Tempo 1 software package performed the fit, 
employing the JPL DE200 solar system ephemeris and the TT(BIPMOl) terrestrial time 
standard of the Bureau International des Poids et Mesures. Orbital kinematics were 
incorporated by means of the theory-independent model of Damour & Deruelle (1986) 
Five Keplerian parameters (orbital period, semi-major axis projected into the line of sight, 
eccentricity, angle of periastron, and time of periastron passage) and one post-Keplerian 
parameter (rate of periastron advance) were necessary to describe the orbit. Also included 
in the timing model were spin period and its time derivative, position and proper motion 
(right ascension and declination and their time derivatives) and a time-varying DM (see 
§3.2). Besides these astrophysical quantities, the fit also allowed for arbitrary time offsets 
between the TOA sets to account for possible alignment discrepancies between standard 
templates as well as for differing signal delays through the various observing hardware. 
Table 2 lists the best-fit timing parameters. 

Included in the table are upper limits on parameters that were not detected. These 
limits were found by allowing the extra parameters to vary one at a time in the fit. In 

addition to the items in the table, we carried out an extensive search for Shapiro delay; see 
§4.4. 


1 http://pulsar.princeton.edu/tempo 
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3.1. TO A Uncertainties 

Preliminary uncertainties of the TOAs were calculated in the profile cross-correlations 
and, for TOAs made by averaging over many short observations, from the scatter in the 
measurements made from the individual short observations. The x 2 of the best fit using 
these uncertainties was high, X 2 /^ = 3.4, where * = 740 is the number of degrees of freedom. 
Most likely this is because the TOA uncertainties had been systematically underestimated, 
a common problem of unknown origin in high-precision pulsar timing. To compensate, we 
added a fixed amount of error in quadrature to the statistical uncertainties of TOAs in each 
data set. The amounts were chosen so that ^jv ~ 1 for each data set in the final fit. 

3.2. Dispersion Measure Variations 

Due to dispersion within the ISM, a radio pulse is delayed in reaching Earth by a 
number of seconds equal to DM/(2.41 x KT 4 / 2 ), where / is the observing frequency in 
MHz and the dispersion measure, DM, is the column density of free electrons integrated 
along the line of sight in units of pc cm 3 , 



in which L is the distance to the pulsar. For many pulsars, the DM can accurately be 
characterized as a single number that holds steady over years of observation. This is not 
true for PSR J0621+1002. Figure la shows the residual pulse arrival times after removing a 
model with a fixed DM. Temporal variations in the DM are visible as secular trends in the 
residuals. As expected for an effect that scales as l// 2 , the lowest frequencies in the figure 
have the largest residuals. We find the pulse profile to be stable across the duration of the 
observations, so that no part of the trends in Figure la results from intrinsic changes in the 
pulsar emission pattern. 
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We incorporated DM variations into the timing model using an 18-term polynomial 
spanning the entire data set, 

17 

DM(<) = DM 0 + - t 0 y, (2 

1=1 

where t 0 is the epoch of the parameter fit and the constant term DM 0 is the value quoted 
for DM in Table 2. The polynomial coefficients were simultaneously fit with all other 
parameters in the global timing solution. The DM 0 term was derived, in effect, from the 
575 and 800 MHz data sets from Green Bank, which were collected with the same data 
acquisition system and timed using the same standard template. The remaining terms 
depended on all the data sets. 

Figure lb shows residual arrival times after subtracting the polynomial DM. These 
residuals are consistent with Gaussian noise except perhaps for an upward rise in the Jodrell 
Bank points in early 2001, a period over which all TOAs belong to a single frequency 
(1400 MHz), so that the DM at those epochs is poorly constrained. These same TOAs have 
some of the largest uncertainties in the data set, however, and so the unmodeled trend in 
them has a negligible effect on the overall timing parameters. 


3.3. Pulsar Astrometry and Spin-down Behavior 

Using the Taylor & Cordes (1993) model of the interstellar free-electron distribution, 
we estimate the pulsar distance to be d — 1.9 kpc. A newer model (J. M. Cordes 2002, 
private communication) puts the pulsar at a distance of d = 1.35 kpc. The spread of 
distance estimates, combined with the total proper motion, // = 3.5 mas yr -1 , gives a 
relative Sun-pulsar transverse velocity in the range V = //d = 23-32km/s. The pulsar is 
at a small distance from the Galactic plane, z = |dsin6| < 0.07 kpc, where b = -2?0 is the 
pulsar’s Galactic latitude. 
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Camilo et al. (2001) have noted that IMBPs tend to have a Galactic scaleheight smaller 
than that of the LMBPs by a factor of 2-4, and that this could be due to their possessing a 
space velocity V smaller than that of the LMBPs by a factor < 2. The low space velocity 
of PSR J0621+1002 supports this notion. With this number there are now 4 IMBPs with 
measured proper motions, for which the space velocities are all approximately 40 ± 10 km s’ 1 
(Toscano et al. 1999; Kramer et al. 1999). These velocities compare to V » 100 km s" 1 for 
a larger population of millisecond pulsars composed largely of LMBPs and isolated pulsars 
(Nice & Taylor 1995; Cordes & Chernoff 1997; Lyne et al. 1998; Toscano et al. 1999). A 
partial explanation for this discrepancy resides in the different evolutionary histories of the 
two binary classes: LMBP progenitors are 1 + 1.3M 0 binaries, while IMBPs may evolve 
from 4 + 1.3 M 0 systems. For identical center-of-mass impulses following the supernova 
explosion, the pulsars in the latter systems will acquire smaller space velocities. 

The observed time derivative of the pulsar’s spin period, P obs , is biased away from 
its intrinsic value, P int , as a result of Doppler accelerations. Following Damour & Taylor 
(1991), we find slight biases due to differential rotation in the plane of the Galaxy and due 
to proper motion. A potential third source of bias, from acceleration perpendicular to the 
plane of the Galaxy, is negligible because of the small distance from PSR J0621+1002 to 
the plane. With the bias subtracted off, we estimate P int = (4.3-4.4) x 10 20 , about 10% 
smaller than P 0 bs- 

Under conventional assumptions about the pulsar spin-down mechanism, the period 
and intrinsic period derivative yield a characteristic age of r = 1.1 x 10 10 yr and a surface 
magnetic field strength of Bq = 1.2 x 10 9 G. 
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4. Pulsar and Companion Masses 

Our goal is to determine the pulsar mass, m j, and the companion mass, m 2 . The 
allowed values of the masses are constrained by the Keplerian orbital elements, the nature 
of the companion star, the apsidal motion of the binary system, and the lack of detectable 
Shapiro delay. In this section, we discuss each of these factors in turn and display the 
resulting constraints on the masses in Figure 2a. The related restrictions on m 2 and orbital 
inclination angle i are shown in Figure 2b. 


4.1. Keplerian Orbit 


The masses and the inclination angle are related through the Keplerian mass function, 


fi(m 1 ,m 2 ,i) 


(m 2 sinz) 3 3 / 27 r \ 2 / 1 \ 

{m l +m 2 y~ x \t~) 


= 0.027026841 ± 0.000000004 M 0 , 


( 3 ) 


where i is the a priori unknown orbital inclination angle; P b is the orbital period; 
x = (aj sin i)/c is the projected semi-major axis of the pulsar measured in light seconds, 
with Oi the semi-major axis and c the speed of light; and T 0 = GM e /c 3 = 4.925 x 10' 6 s, 
with G Newton’s gravitational constant. Since sin* < 1, this can be rewritten to give an 
upper limit on m x in terms of m 2 , 


mi < m^ 2 f 1 1/2 - m 2 . 

This constraint is shown as the lower shaded region in Figure 2a. 


( 4 ) 


4.2. Upper Limit on Companion Mass 

The companion to the pulsar must be either a main sequence star, a neutron star, a 
black hole, or a white dwarf. The first possibility can be ruled out, as optical observations 
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using the Hubble Space Telescope find no evidence that the secondary is a main sequence 
star (Camilo et al., to be published elsewhere). The second possibility can also be 
eliminated, since it is improbable that a double neutron star binary could have survived 
two supernova explosions and yet retained the small eccentricity of the PSR J0621+1002 
system (Portegies Zwart & Yungelson 1998). A pulsar-black hole binary is also unlikely to 
have such a circular orbit (Lipunov et al. 1994). The companion must therefore be a white 
dwarf, and as such, its mass, m 2 , must be lower than the Chandrasekhar limit, 1.4 M©. 
This constraint is illustrated in Figure 2 as the upper shaded regions. 


4.3. Relativistic Periastron Advance 


The measured rate of periastron advance, u, provides another relation between m\ 
and m 2 . For the PSR J0621+1002 system, we assume that nonrelativistic contributions to 
apsidal motion are negligibly small (see §4.6). The general relativistic interpretation of w 
then yields the combined mass of the stars, 


M = mi + m 2 



5/2 f(l -e 2 )iJ 


3/2 



2.81 ±0.30 M 0 . 


( 5 ) 


This constrains m i and m 2 to lie within the strips indicated by dashed lines in Figure 2a. 


4.4. Shapiro Delay 

According to general relativity, a pulse is delayed as it propagates through the 
gravitational potential well of the secondary. For a pulsar in a circular orbit, this so-called 

Shapiro delay is 


A t s = -2m 2 T Q ln[l - sinisin^ - </>o)]> 


( 6 ) 
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where 0 is the orbital phase in radians and 0 O is the phase of the ascending node. For 
small inclination angles, the variation of A t s over an orbit is nearly sinusoidal, and so it is 
indistinguishable from a minor increase in the projected orbital size, x. For edge-on orbits, 
with i « 90°, the variation becomes strongly peaked at 0 - « tt/ 2, when the pulsar is 

behind the companion. This breaks the covariance with the Keplerian parameters, allowing 
measurement of Shapiro delay, and hence of 77^2 and i . 

Shapiro delay is not detected in the PSR J0621+1002 timing data. However, the 
magnitude of the effect is expected to be around 15 //s, which would make Shapiro delay 
easily detectable if the inclination angle were large. The fact that the delay is not observed 
therefore implies that the orbit is tilted substantially away from an edge-on orientation. 

To explore the statistical limits that both the detection of periastron advance and the 
non-detection of Shapiro delay place on i and m 2 (and hence on m,), we analyzed the 
timing data over a grid of values in the ranges 0 < cos i < 1 and 0 <m 2 < 2.0 M 0 . For each 
combination of m 2 and cosi, we calculated the Shapiro delay parameters and the rate of 
periastron advance according to general relativity. We then performed a timing fit holding 
those quantities fixed while allowing all other parameters to vary. We recorded the resulting 
value of x 2 and its difference from the global minimum on the grid. Small departures in y 2 
from the minimum signified the most likely configurations of m 2 and cos i. 

The solid contours in Figure 2a show the regions in which acceptable timing solutions 
were found. The interpretation of the contours is straightforward: within the area allowed 
by the precession measurement, the non-detection of Shapiro delay excludes solutions with 
edge-on orbits, and so the strip is truncated at low values of cos i. Figure 2b shows that 
the inclination angle is constrained with a high degree of confidence to be less than 50°. 

Following the statistically rigorous procedure described in Appendix A, we converted 
the x 2 differences to probabilities and derived probability distribution functions (PDFs) 
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for mi and m 2 . The analysis is restricted to timing solutions for which m 2 < 1.4 M 0 , in 
accordance with the discussion in §4.2. The PDFs of the pulsar and companion masses 
appear in Figure 3. They yield the maximum likelihood estimates 

mi = 1.69152 M 0 (68% confidence) ( 7 ) 

for the pulsar and 

m 2 = 0.98152 M © (68% confidence) (8) 

for the companion. Note that the sum of the maximum likelihood estimates, 2.67 M©, is less 
than the total system mass derived in equation 5 from periastron advance alone, 2.81 M 0 . 
This is primarily a consequence of the upper limit on the companion mass, m 2 < 1.4 M Q , 
which preferentially excludes solutions with high total mass, as can be seen in Figure 2a. 


4.5. Interpretation of the Masses 

A growing body of evidence finds that neutron stars in white dwarf-pulsar systems 
are not much more massive than those in double neutron star binaries, even though the 
secondaries must lose several tenths of a solar mass as they evolve toward white dwarfs. In 
particular, Thorsett & Chakrabarty (1999) find that the masses of neutron stars orbiting 
either white dwarfs or other neutron stars are consistent with a remarkably narrow Gaussian 
distribution, m = 1.35 ± 0.04 M©. While our measurement of the mass of PSR J0621+1002 
is in statistical agreement with that result, the maximum likelihood value is suggestively 
high, allowing the possibility that a substantial amount of mass was transferred onto the 

neutron star. 

Our estimate of the companion mass implies that the star is probably an ONeMg 
white dwarf and ranks it among the heaviest known white dwarfs in orbit around a pulsar. 
As van Kerkwijk k Kulkarni (1999) have pointed out, there are few massive white dwarfs, 
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m 2 >1M 0 , known to have evolved from a single massive progenitor star. They argue that 
the companion to PSR B2303+46, a young pulsar in an eccentric orbit, is such a white 
dwarf. Our timing of PSR J0621 + 1002 shows that its companion is similarly likely to 
have descended from a massive star, even though the histories of mass loss and accretion 
in the eccentric PSR B2303+46 binary, with no spin-up of the pulsar, and the circular 
PSR J0621 + 1002 binary, with significant spin-up, must have been very different. 

Given the mass estimates, it seems likely that the PSR J0621+1002 system formed 
through a common envelope and spiral-in phase (Taam, King, & Ritter 2000; Tauris, van 
den Heuvel, & Savonije 2000). In this scenario, the companion originally had a mass of 
5 — 7M 0 , and the pulsar, initially in a wide orbit with a binary period of a few hundred days, 
spiraled in to its current orbit of P b = 8.3 under a drag force arising from its motion through 
the envelope. The same formation mechanism has been put forward for PSR J1454-5846, 
an object quite similar to PSR J0621+1002, with a 12.4-day orbital period, a 45.2-ms spin 
period, an eccentricity of 0.002, and a companion mass of ~ 1.1 M 0 (Camilo et al. 2001) 

4.6. Classical Periastron Advance 

The above analysis assumed that the observed apsidal motion can be entirely attributed 
to relativity. In principle, however, apsidal motion could also be caused by distortions of 
the secondary star. Smarr & Blandford (1976) considered this possibility in the context 
of a potential white dwarf companion to the Hulse-Taylor binary pulsar, PSR B1913+16 
Their analysis can also be applied to PSR J0621+1002. They found tidal deformation of 
the companion to contribute negligibly to u for PSR B1913+16. Since the apsidal advance 
per binary period due to tidal deformation scales as a~ 5 , where a is the major axis, this 
effect can also be neglected in the wider PSR J0621+1002 system. 
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A potentially more important effect is rotational deformation, which becomes significant 
if the secondary is spinning rapidly. The precession rate due to rotation is 

throt = nQ (l - | sin 2 9 + cot i sin 6 cos 0 cos <F 0 ) ( 9 ) 

(Wex 1998), in which n = 2n/P b and 

3 fez-Rffi 2 (10) 

a*Gm 2 ( 1-e 2 ) 2 ’ 

where k 2 , # 2 , ^ 2 , and m 2 are the structure constant, radius, angular velocity, and mass 
of the secondary; 9 is the angle between the angular momentum vector of the secondary 
and the angular momentum vector of the orbit; and $0 is the longitude of the ascending 
node in a reference frame defined by the total angular momentum vector (see Fig. 9 of Wex 
1998). Neither 9 nor $0 is known, so we must allow for all possible values in the ranges 
0 < 9 < 7 T and 0 < $ 0 < 2 tt. Substituting PSR J0621+1002’s Keplerian orbital parameters, 

the formula can be written in the notation of Smarr & Blandford (1976). 

( A f \ 2/3 

m'+m ) a6 l 1 ~ ^ Sin2 6 + cot * sin 6cOS dcos ®°) ’ 

where a = 2k 2 H^Q 2 2 /(ZGm 2 ) and a 6 = a/(10 6 km 2 ). To gauge the largest value that 
<j rot could attain, we will use limits on the masses and the orbital inclination angle from 
the relativistic analysis, recognizing that they would need to be modified should u rot be 
found significant. Our observations constrain 20° < i < 50°, so that 1.2 < coti < 2.7 
(see Figure 2b), for which the maximum value of the geometric factor in equation 11 is 
1.7, attained at 9 — 40° and <h 0 = 0. Analyzing models of rotating white dwarfs, Smarr 
& Blandford (1976) found that a 6 < 15. Combining these restrictions with our measured 
value of total system mass, the precession due to rotation for PSR J0621+1002 can be as 
high as cj rot < 0?0021 yr~\ about 20% of our observed value. 

There is reason to believe that the classical contribution to the observed u is smaller 
than this. In most cases, a rotational deformation will induce a change in the projected 
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semi-major axis of the orbit, x = a x sin i/c. Wex (1998) shows the rate of change to be 


from which 


^rot — 


vi 

~ —nQ cot i sin 9 cos 9 sin 4> 0 , 
x 


x / 1 - | sin 2 9 + cot i sin 6 cos 0 cos $ 0 


) 


( 12 ) 


(13) 


x \ cot i sin 9 cos 9 sin <h 0 

Thus, our observed upper limit of \x/x\ < 1.2 x 10~ 15 implies that w rot is no more than 
1.2 x 10 -15 rads -1 = T.l x 10“ 5 6 yr _1 times a geometric factor. Unfortunately, certain 
special combinations of i, 0, and $o will make the geometric factor large. For this reason, 
we cannot definitively exclude the possibility that rotational precession contributes to the 
observed value of Co. For most cases, however, the geometric factor will be of order unity, 
and so the upper limit on ch ro t will be substantially smaller than our observed value of 
Co. Because of this, and because there is no reason to expect the secondary to be rapidly 
rotating, we have chosen to ignore w rot in our analysis of the pulsar and companion masses. 


5. Density Irregularities in the ISM 
5.1. Temporal DM Variations 

Figure 4a shows the DM of PSR J0621+1002 calculated on individual days on which 
data were collected at Arecibo at both 430 and 1410 MHz. A DM drift as steep as 0.013 
pc cm 3 yr 1 can be seen in the plot. If not properly modeled, such a gradient would 
shift the TOAs by up to 7 fis at 430 MHz over the 8-day pulsar orbit, a systematic effect 
significantly larger than the measurement uncertainties of the Arecibo TOAs. The DM 
gradient is among the largest ever detected in a pulsar outside a nebula (Backer et al. 1993) 

Variations in the DM can be used to investigate inhomogeneities in the interstellar 
electron density. For this purpose, we modeled the DM as a series of step functions in 
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time, as shown in Figure 4b. The width of the step intervals, 100 days, was a compromise 
between the goal of sampling DM as frequently as possible and the need for each segment 
to span enough multi-frequency data for the DM to be reliably calculated. 

In Figure 4b, there appear to be discontinuities in the DM at 1997.4 and 1999.0. With 
a careful check of the data around these dates, we confirmed that the DM did indeed 
change by the amounts shown within the 100-day time resolution of the figure, and that 
the increase and decrease are not artifacts of binning the DM or of joining together data 
from different receivers and telescopes. However, the data do not allow us to distinguish 
between nearly instantaneous jumps in DM, as seen in the Crab pulsar signal (Backer, 
Wong, & Valanju 2000), and slower changes on a scale of 100 days. Given the rapid but 
apparently smooth variations in DM seen after 2000, we suspect the DM to be strongly but 
continuously varying at the earlier epochs as well. 


5.2. Structure Function Analysis 

The DM variations in Figure 4 result from the passage across the line of sight of density 
irregularities in the ionized ISM. The spatial structure of the irregularities can be discerned 
using the two-point structure function of the DM, 

D d „(t) = <[DM(! + r) - DM(f)] 2 ) . (14) 

where r is the time lag between DM measurements and where the angular brackets denote 
ensemble averaging (e.g., Cordes et al. 1990). In a simple one-dimensional model where 
the line of sight cuts at the relative Sun-pulsar transverse velocity V across a pattern of 
ISM irregularities that are “frozen in” a thin screen, and where the thin screen is midway 
between the Sun and the pulsar, the time lag is related to the spatial size l of the density 
inhomogeneities through l = Vt/2. 



The structure function ■Ddm(t') can be approximated from the DM values at times t,- 
through the unbiased estimator 

1 n(t} 

^Dm(t) = N\t) + r ) - DM (^)] 2 - a DM, (15) 

where N(t) is the number of DM pairs that enter into the summation at a given lag r and 
where <r DM is the mean uncertainty of the DM values. Figure 5 illustrates values of D DM ( r ) 
calculated this way for PSR J0621+1002. The error bars were computed by assuming 
Gaussian statistics for the fitted DM values and formally propagating their covariances 
through equation 15. We restrict r to values for which N(t ) > 10, thereby extending the 
lag from a minimum of 100 days to a maximum of 1300 days, about half the length of our 
full data span. 


5.3. Not a Simple Power Law 

The structure function in Figure 5 presents a puzzle because it is not the simple power 
law that is predicted by standard theories of the ISM and that is seen in the direction 
of other pulsars. In the standard picture, a power law arises from the expectation that 
turbulence spreads energy from longer to shorter length scales. Accordingly, the spectrum of 
perturbations in the electron density is modeled in terms of the spatial frequency, q = 2n/l, 
via 

P(q) oc q ~ 0 (16) 

(Rickett 1990). The relation is hypothesized to hold over some range of wave numbers, 

Qi < Q < q a \ below some “inner scale” q t = 2i x/U and above some “outer scale” q Q = 2i x/l 0l 
various damping mechanisms are expected to cut off the turbulent energy flow. With the 
assumed linear relation between l and t, the structure function becomes 
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where tq is 9. normalization constant. The index j3 is usually predicted to be near the 
Kolmogorov value for turbulence in neutral gases, 11/3. Studies of DM variations in other 
pulsars have uncovered power laws with indices close to that value (Phillips k Wolszczan 
1992; Kaspi, Taylor, k Ryba 1994; Cognard k Lestrade 1997). 

The phase structure function in Figure 5 is obviously not a simple power law. It is not 
clear how to interpret this. A power law, marked by the dashed line, can be fit to the first 
five points in the plot. The slope of the line yields a spectral index, (3 3.28 i 0.09, that 

is reasonably close to the Kolmogorov value. At lags longer than 500 days, the power law 
fails to hold, although there is a hint of its reemergence at the longest lags. The flattening 
out of the structure function at 500 days is a direct consequence of the two sharp breaks 
in the DM time series (see Fig. 4b). If we take the Sun-pulsar transverse velocity to be 
V = 27kms~ 1 , then 500 days corresponds to a length of l = Vt/2 = 5.5 x 10 13 cm, implying 
structure in the electron density power spectrum at this scale. 


6. Conclusion 

We have found substantial constraints on the masses of PSR J0621+1002 and its 
orbital companion. The pulsar mass is found to be m\ = 1.69^ 3 0 M©- The lower end of 
this uncertainty range is near the canonical pulsar mass of 1.35 M©, but the mass may also 
be several tenths of a solar mass higher, allowing the possibility that a substantial amount 
of material accreted onto the neutron star during the evolution of the system. The mass of 
the secondary, m 2 = 0.98+°; 3 2 M©, makes it one of the heaviest known white dwarfs orbiting 
a pulsar. 

Can the mass measurements be improved by continued timing observations? 
Unfortunately, post-Keplerian effects beyond those considered in this paper, such as orbital 
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period decay, and gravitational redshift and time dilation, will not be detectable in the 
timing data for PSR J0621+1002 in the foreseeable future, so any improvement must come 
about through improved measurements of u> and Shapiro delay. 

The uncertainty in the measurement of the total mass, M = 2.81 ± 0.30 M©, scales 
linearly with the uncertainty of tu, which in turn is inversely proportional to the time 
span over which data are collected. This has a simple explanation: the longer the time 
span of the observations, the more u shifts, and so the easier it is to measure c j. The 
highest precision data used in this work — the annual Arecibo campaigns — were collected 
over two years. A similar campaign carried out several years in the future would shrink the 
uncertainty in M by a factor of a few. 

The pulsar and companion masses were further constrained by Shapiro delay. The 
precision of the Shapiro delay measurement (or limit) has no dependence on data span 
length, so its uncertainty is reduced only as n -1 / 2 , where n is the number of observations. 
At best, a modest improvement could be made with existing radio telescope resources. 
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timing program; we particularly thank D. Backer, A. Lommen, D. Lorimer and K. Xilouris 
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Center, which is operated by Cornell University under a cooperative agreement with the 
National Science Foundation. The National Radio Astronomy Observatory is operated by 
Associated Universities, Inc., under a cooperative agreement with the National Science 
Foundation. Pulsar research at Princeton University is supported by National Science 
Foundation grant 96-18357. FC is supported by SAO grant GO1-2063X. IHS is a Jansky 
Fellow. 
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A. Statistical Analysis of Pulsar and Companion Masses 

The PDFs of the pulsar mass, m u and the companion mass, m 2 , shown in §4.4 were 
calculated by least-squares fitting the TOAs to a timing model across a grid of m 2 and cos z 
values and analyzing the resulting changes in x 2 - A formal derivation of the PDFs proceeds 
as follows. If Xo is the global minimum on the grid, then each value of 

Ax 2 (m 2 , cosz) = x 2 ( m 2 , cos i) - Xo (Al) 

has a x 2 distribution with two degrees of freedom. It therefore maps to a Bayesian likelihood 
function, 

p{{tj} | m 2 , cost) = Q'j e" Ax2/2 , (A2) 

where {ij stands for the data set. Accordingly, the joint posterior probability density of 
m 2 and cosz is 

p(m 2 ,cosz | {*,}) = — - ^ ^° — p(m 2 , cos i) . (A3) 

The Bayesian “evidence”, p({tj}), is determined by normalizing the integral of 
p(m 2 ,cosz |{tj}) over all grid points that are consistent with mi > 0, given the relationship 
among mi, m 2 and cosz in the mass function in equation 3. As in all Bayesian investigations, 
a choice must be made for the prior, p(m 2 ,cosz). We selected the product of a uniform 
distribution on 0 < m 2 < ra 2)max and a uniform distribution on 0 < cosz < 1. For the 
companion mass, the choice embodies our ignorance about the star; all that is known for 
certain is that it is a white dwarf, for which reason we set m 2i max 1*4 M 0 (see §4.2). For 
the inclination angle, a flat distribution in cos z follows from the assumption that the orbital 
angular momentum vector has no preferred direction in space. 

The PDF of m 2 is obtained by marginalizing equation A3 with respect to cosz: 

p(m 2 \{tj}) = f cl{cos i)p(m 2 , cos i\{tj}). (A4) 

Jo 
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Similarly, the PDF of m\ can be expressed as a double integral, 

prn 2, max pi 

p{m^\{tj}) = / dm <2 / d(cosi)p(mi \m 2 , cos i) p(m 2 , cos i |{t,}), 
Jo Jo 

where mi is guaranteed to be consistent with the mass function by setting 
p(m 1 |m 2 , cos?) = 5 


(A5) 


(m 2 sin i) 3 / 2 

™i ~ I ~T7o rn 2 


fi 


1/2 


(A6) 


with S the Dirac function. 
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Table 1: Summary of Observations. 


Observatory 

Dates 

Frequency 

(MHz) 

Bandwidth 

(MHz) 

Number of 
TOAs 

Typical Integration 
Time (min) 

RMS residual* 
(ns) 

Arecibo 

1997.9-2001.5 

430 

5 

103 

29 

2.6 


1999.4-2001.5 

1410 

10 

24 

29 

3.2 

Green Bank 

1995.2-1999.0 

370 

40 

48 

50 

17 


1995.2-1999.5 

575 

40 

49 

40 

13 


1995.2-1999.5 

800 

40 

39 

40 

19 

Jodrell Bank 

1996.2-1997.7 

410 

40 

12 

30 

21 


1995.7-1999.8 

606 

40 

298 

30 

12 


1995.7-1997.8 

1410 

32 

78 

30 

24 


1997.8-1999.3 

1380 

96 

79 

30 

16 


1999.7-2001.2 

1396 

64 

50 

30 

19 


“Values for Arecibo and Green Bank include the effect of averaging TOAs calculated from shorter integration times. 
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Table 2: Pulse Timing Parameters of PSR J0621 + 1002 a 

Measured Parameters 


Right ascension, a (J2000) 06 h 21 m 22? 11108(3) 

Declination, 8 (J2000) + 10° 02' 38'/ 74 1(2) 

Proper motion in R.A., p a = a cos 8 (mas yr" 1 ). , . 3.5(3) 

Proper motion in Dec., =8 (mas yr -1 ) -0.3(9) 

Pulse period, P (ms) 28.853860730049(1) 

Period derivative, P 0 ^ s (10" 2 °) 4.732(2) 

Epoch (MJD [TDB]) 50944.0 

Orbital period, P b (days) 8.3186813(4) 

Projected semi-major axis, x (lt-s) 12.0320744(4) 

Eccentricity, e 0.00245744(5) 

Epoch of periastron^To (MJD [TDB]) 50944.75683(4) 

Longitude of periastron^u; (deg) 188.816(2) 

Periastron rate of change, uj (deg yr" 1 ) 0.0116(8) 

Dispersion measure? DM (pc cm" 3 ) 36.6010(6) 

Measured Upper Limits 

Parallax (mas) < 2.7 

Pulse period second derivative, P (s" 1 ) <4x 10" 3 ° 

Orbital period rate of change, P b < 5 x 10" 12 

Orbital axis rate of change, x < 1.5 x 10" 14 

Derived Parameters 

Mass function, /i (M©) 0.027026841(4) 

Total mass, M (M©) 2.81 ± 0.30 

Pulsar mass, mi (M©) 169 lo30 

Companion mass, m 2 (M©) 0.981^2 

Characteristic age (yr) 1.1 x IO 10 

Surface magnetic field strength (Gauss) 1.2 x 10 9 

Total proper motion, /i (mas yr" 1 ) 3.5(3) 


“Values in parentheses are la uncertainties (68% confidence) in the last digit quoted. 

b u and To are highly covariant. Observers should use u; = 188?815781 and T 0 = 50944.756830176. 

c The DM varies. The value here is the constant term in an 18-term polynomial expansion (see §3.2). 
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E1 


E1 NOTE TO EDITOR: Please place Figures 2a and 2b side by side instead of on top of 


one another. 
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(a) Fixed DM 


(b) Variable DM 
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Fig. 1. — Residual pulse arrival times for fits with (a) the DM held constant and (b) the DM 
modeled as an 18-term polynomial. Note the difference in the vertical scales of the two plots 
and, in particular, the scale of the bottom Arecibo panel, emphasizing the high precision of 
the TOAs from that observatory. 
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Fig. 2. — Values of pulsar mass, mi, companion mass, m 2 , and orbital inclination angle, 
cosi, permitted by the Keplerian orbital elements, the relativistic timing model and the 
assumption that the companion is a white dwarf, (a) Allowed values of mi and m 2 . Dotted 
lines show selected values of constant i. (b) Allowed values of cos i and m 2 . In both plots, 
the true values of the parameters must lie within the A% 2 contours and outside the gray 


regions. 
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Fig. 3.— Probability distribution functions of (a) the pulsar mass and (b) the companion 
mass. The solid vertical lines mark the centers of the peak histogram bins. Their locations 
are quoted in the text as the maximum likelihood values of the masses. The dashed lines are 
error bars around those values. They represent the shortest widths along the mass axes that 
both enclose the peaks and contain 68% of the area under the curves. For the companion 
mass, the sharp falloff in the histogram reflects the assumption that m 2 < 1.40 M 0 . 
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Fig. 4. — Dispersion measure as function of time, with ADM = (DM — 36.6010). Vertical 
error bars are la uncertainties in the variations, (a) Measurements within a single day 
derived from dual-frequency observations at Arecibo. (b) Measurements derived from all 
TO As averaged over intervals of 100 days, as indicated by horizontal bars. 






- 32 - 


h 


Q 


<Q 

i i 

tuO 

O 



r (days) 


Fig. 5. — Structure function of DM variations of PSR J0621+1002 (Fig. 4). The points 
clearly do not all follow a single power law. As indicated by the dashed line, a power law 
was least-squares fitted to the first five points according to equation 17, yielding the spectral 
index, j3. 




